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_ We consider a stochastic volatility model with Levy jumps for a log-return process Z = (Zt)t>o of the form 

Pj^ ■ Z — U + X , where U = (Ut)t>o is a classical stochastic volatility process and X — (X t )t>o is an independent 

Levy process with absolutely continuous Levy measure v. Small-time expansions, of arbitrary polynomial order, in 
time-t, are obtained for the tails P (Z t > z), z > 0, and for the call-option prices E (e z+Zt - l)^, z / 0, assuming 
smoothness conditions on the density of v away from the origin and a small-time large deviation principle on U. 
Our approach allows for a unified treatment of general payoff functions of the form (p(x)l x > z for smooth functions 
J"V^ tp and 2 > 0. As a consequence of our tail expansions, the polynomial expansions in t of the transition densities ft 

0^ ■ are also obtained under mild conditions. 
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1 Introduction 



It is generally recognized that accurate modeling of the option market and asset prices requires a mixture of a con- 
tinuous diffusive component and a jump component. For instance, based on high-frequency statistical methods for 
Ito scmimartingales, several empirical studies have statistically rejected the null hypothesis of either a purely-jump 
or a purely-continuous model (see, e.g., [T], [J, [5], [B], [35] ). Similarly, based on partially heuristic arguments, [TU] 
characterizes the small-time behavior of at-the-money (ATM) and out-of-the-money (OTM) call option prices, and 
based on these results, then argues that both, a continuous and a jump component, are needed to explain the behavior 
of the market implied volatilities for S&P500 index options. Another empirical work in the same vein is [30] . where 
a small-time small-log-moneyness approximation for the implied volatility surface was studied in the case of a local 
jump-diffusion model with finite jump activity. Based on S&P500 option market data, [30j also concludes that jumps 
are significant in the risk- neutral price dynamics of the underlying asset. 

Historically, local volatility models (and more recently stochastic volatility models) were the models of choice to 
replicate the skewness of the market implied volatilities at a given time (see [19] and [21] for more details). However, it 
is a well-known empirical fact that implied volatility skewness is more accentuated as the expiration time approaches. 
Such a phenomenon is hard to reproduce within the purely-continuous framework unless the "volatility of volatility" is 
forced to take very high values. Furthermore, as nicely explained in (Chapter 1), the very existence of a market for 
short-term options is evidence that option market participants operate under the assumption that a jump component 
is present. 

In recent years the literature of small-time asymptotics for vanilla option prices of jump-diffusion models has grown 
significantly with strong emphasis to consider either a purely-continuous model or a purely Levy model. In case of 
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stochastic volatility models and local volatility models, we can mention, among others, [7], jS], [12], [17], [H], [20], [23], 
[33] . For Levy processes, [35] and [3j5] show that OTM option prices are generalljQ asymptotically equivalent to the 
time-to-maturity t as t — > 0. In turn, such a behavior implies that the implied volatilities of a Levy model explodes 
as r — > 0. The exact first order asymptotic behavior of the implied volatility close to maturity was independently 
obtained in j!4) and [39] . while the former paper also gives the second order asymptotic behavior. There are very 
few articles that consider simultaneously stochastic volatility and jumps in the model. One such work is [10] which 
obtains, partially via heuristic arguments, the first order asymptotic behavior of an Ito semimartingale with jumps. 
Concretely, ATM option prices of pure-jump models of bounded variation decrease at the rate O(t), while they are 
just 0(i/r) under the presence of a Brownian component. By considering a stable pure-jump component, [10] also 
shows that, in general, the rate could be O(r^), for some j3 € (0, 1). For OTM options, it is also argued that the first 
order behavior is 0(e~ c / T ) in the purely-continuous case, while it becomes 0(t) under the presence of jumps. Very 
recently, [31] formally shows that the leading term of ATM option prices is of order \ T for a relatively general class of 
purely-continuous Ito models, while for a more general type of Ito processes with a- stable- like small jumps, the leading 
term is 0{T X l a ) (see also [14] Proposition 4.2] and [39] Proposition 5] for related results in exponential Levy models). 
Fractional expansions are also obtained for the distribution functions of some Levy processes in |29| . 

In this article, we consider a jump diffusion model by combining a stochastic volatility model with a pure-jump 
Levy process of infinite jump activity. More precisely, we consider a market consisting of a risk-free asset with constant 
interest rate r > and a risky asset with price process (St)t>o defined on a complete probability space (f2,J-", Q) 
equipped with a filtration (.Ft)t>o satisfying the usual conditions. We assume that, under Q, the log-return process 

Z t :=\o^— \ (1.1) 

DO 

follows the jump diffusion model 

Zt = U t + X t , dUt = fx(Y t )dt + a(Y t )dW t (1) , (1.2) 
dY t = a(Y t )dt + j(Y t )dW t {2) (1.3) 

with U = X = 0, Y = y e K. Here, and are Wiener processes adapted to (J 7 *), X is an independent 

(J r t)-adapted pure-jump Levy process with triplet (b,0,v), and er, 7, /1 and a are deterministic functions such that 
(|1.2[) - (|1.3[) admits a solution. We also assume that f 1 e z v(dz) < 00 and 

e z -1- zl\ z \< 1 )v{dz), and fi(y) = -\<J 2 {y). (1.4) 

In particular, note that Q is assumed to be a martingale measure (i.e. (e~ rt St)t>o is a Q-martingale). The model 
(|1.2p - (|1.3j) is appealing in practice since it incorporates jumps in the asset price process as well as volatility clustering 
and leverage effects. The process (Yt)t>o is the underlying volatility factor driving the stochastic volatility of the 
process. 

For z ^ and t > 0, let 

G t (z):=E(e z+z *-l) + , (1.5) 

where E denotes, from now on, the expectation under Q. We will show that, under mild conditions, the following 
small-time expansions for Gt hold true: 



G t (z) = 5>-(z)- + 0(t" +1 ), (1.6) 



3=0 J 



for each n > and certain functions by Note that the time-i price of a European call option with strike K, which is 
not at-the-money, when the spot stock price and the underlying volatility factor have respective values s and yo, can 
then be expressed as 

C{t,s) :=e- r ( T -* ) E((S , T - K) \ F t ) = Ke TT G T (lns — hiK), (1.7) 



lr That is, except for some pathological cases (see 1361 for examples) 
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where r = T — t and s ^ K. Hence, the small-time behavior of (|1.5[) leads to close-to-expiry approximations for the 
price of an arbitrary not-at-the- money call option as a polynomial expansion in time-to-maturity. From note 
also that the expression 

S e- z G t (z) = e- rt E (S t - S e rt - Z ) + (1.8) 

can be interpreted as the time-i call option price with log-moneynes^| n := —z. 

Small-time option price asymptotics for the model (|1.2[) - (|1.3[) were also considered in [30], but only for finite-jump 
activity Levy processes X. Another important difference is that, here, we focus on small-time asymptotics for fixed 
log-moneyness z / 0, while |30j considers approximations where z is simultaneously made to converge to as t — > 
(small-time and small- log-moneyness asymptotics). Let us also remark that [30] assumes throughout (and without 
proof) that the implied volatility surface satisfies an expansion in powers of z and t, which, in our opinion, is a rather 
strong assumption. Another related work is |31j . where the first order small-time small-log-moneyness asymptotics 
is considered for a class of Ito semimartingales with non-zero continuous component (Theorem 3.1 therein). For the 
CGMY and related models, higher order approximations for ATM option prices are obtained in [T5] , 

Our method of proof is built on a type of iterated Dynkin formula of the form 

™ xk j.n+1 A 

Eg(X t ) = g(0) + V T .L k g$) + — r / (1 - a)"E (L n+1 g{X at )) da, (1.9) 
fe=1 fc! n. J 

where g is a sufficiently smooth function and L is the infinitesimal generator of the Levy process X given by 

2 p 

Lg{x) := ys"(aO + W(x) + J {g(z + x)- g(x) - zg'(x)l {[z{ < 1} ) v(dz), (1.10) 

for g £ Cl and a Levy triplet (b, a 2 , v) (see Section [2TT1 below for terminology). The main complication with option call 
prices arises from the lack of smoothness of the payoff function g z (x) = (e z+x — 1)+. In order to "regularize" the payoff 
function g, we follow a two step procedure. First, we decompose the Levy process into a compound Poisson process 
with a smooth jump density vanishing in a neighborhood of the origin and an independent Levy process with small 
jumps. Then, we condition Eg(X t ) on the number of jumps of the compound Poisson component of X and apply the 
Dynkin's formula on each of the resulting terms. Contrary to the approaches in [14| and [39] . where the special form 
of the payoff function g z (x) = {e z+x — 1)+ plays a key role, our approach can handle more general payoff functions of 
the form 

gz(x) = ip(x)i {x > z} , (l.n) 

for a smooth function ip and positive z. By considering the particular case ip(x) = 1, we generalize the distribution 
expansions in |13] to our jump-diffusion setting. Let us emphasize that the process Z in (jl.2|) is not truly a Markov 
model but rather a hidden Markov model. This fact causes some technical subtleties that require a careful analysis of 
the iterated infinitesimal generator of the bivariate Markov process {(U t , Y t )} t >o- 

As an equally relevant second contribution of our paper, we also obtain polynomial expansions for the transition 
densities ft of the Levy process, under conditions involving only the Levy density of X . This is an important improve- 
ment to our former results in [13j . where a uniform boundedness condition on all the derivatives of ft away from the 
origin was imposed. Expansions for the transition densities of local volatility models (with possibly jumps but only of 
finite activity) have appeared before in the literature (e.g. see [2], [3], [40]). Unlike our approach, the general idea in 
these papers consists of first proposing the general form of the expansion, and then choosing the coefficients so that 
either the backward or forward Kolmogorov equation is satisfied. The resulting coefficients typically involve iterated 
infinitesimal generators as in our expansions, even though our approximations are uniform away from the origin. 

The paper is organized as follows. Section 2 contains some preliminary results on Levy processes, which will be 
needed throughout the paper. Section 3 establishes the small-time expansions, of arbitrary polynomial order in t, for 
both the tail distributions F(Z t > z), z > 0, and the call-option price function Gt(z), z ^ 0. This section also justifies 
the validity of our results for payoff functions of the form (jl.llj) . Section 4 illustrates our approach by presenting 
the first few terms of those expansions. Interestingly enough, the first two coefficients of the expansion of the general 
model coincide with the first two coefficients of an exponential Levy model. Section 5 obtains the asymptotic behavior 
of the corresponding implied volatility. Section 6 gives a small-time expansion for the transition density of a general 
Levy process under rather mild conditions. The proofs of our main results are deferred to Appendices. 

2 As usual, log-moneyness is denned as the logarithm of the ratio between the strike K and the forward price Soe rt . 
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2 Background and preliminary results 



2.1 Notation 



Throughout this paper, C n (or C n (R)), n > 0, is the class of real valued functions, defined on K, which have 
continuous derivatives of order < k < n, while C C" corresponds to those functions having bounded derivatives. 
In a similar fashion, C°° (or C°° (R)) is the class of real valued functions, defined on R, which have continuous derivatives 
of any order k > 0, while C£°(R) C C°° arc again the ones having bounded derivatives. Sometimes R will be replaced 
by Ro := R \ {0} or R fc when the functions are defined on these spaces. 

Throughout this section, X denotes a Levy process with triplet (6, er 2 , v) defined on (Q, J 7 , (J r t)t>o, Q)- Let us write 
X in terms of its Levy-Ito decomposition: 

X t = bt + aW t + [ [ zn(ds,dz)+ f f zp(ds,dz), 

Jo J\z\>l JO J\z\<l 

where W is a Wiener process and \i is an independent Poisson measure on R + x R\{0} with mean measure dtv(dz) 
and compensator /2. For each e > 0, let c £ € C°° be a symmetric truncation function such that lr_ e /2 e /2] ( z ) < c s( z ) ^ 
lr_ e>e ](z) and let c e := 1 — c s . Next, for < e < 1, consider two independent Levy processes, denoted by X £ and X £ , 
with respective Levy triplets (0, 0, c e (z)v{dz)) and (b £ , a 2 , c e {z)v{dz)), where 



zc e (z)z^((i2;). 



z <i 



Note that (Xt)t>o has the same law as (X| + Xf)t>o and, therefore, without loss of generality, we assume hereafter 
that X = X s + X e . Note also that X e is a compound Poisson process with intensity of jumps X £ := J c e {z)v{dz) and 
jumps distribution c e (z)u(dz)/X £ . Throughout, (£f)i>i stands for the jumps of the process X s . The process X e has 
infinitesimal generator L e given by 

2 

L s g(y) = b e g'{y) + ^-g"(y) + l e g(y), (2.1) 

for g G C 2 , where 

2e5(2/) := / {.9(2/ + z) - - zg' {y)l\ z \< x } c e {z)v{dz) 
g"(y + Pw)(l - P)d/3w 2 c e {w)v{dw). 



/Ro JO 

The following tail estimate for X £ is also used in the sequel: 

P(|X t e | > z) < t az exp(az lnz )exp(az - azhaz), (2.2) 

where a £ (O,^ 1 ), and t, z > satisfy t < z/z for some zq depending only on a (see [351 Section 2.6], [371 Lemma 
3.2] and [T3J Remark 3.1] for proofs and extensions). 

Throughout the paper, we also make the following standing assumptions: 

v{dz) = s{z)dz, seC*°°(R\{0}) and 7 M := sup |s (fc) (z)| < 00, for all 8 > 0, (2.3) 

|z|>a 

e c ^v{dz) < 00, for some c > 2. (2.4) 

\z\>l 

Finally, the following notation are also in use: 

s £ :=c e s, s £ :=(l-c £ )s, L°g = g, L k+1 g = L{L k g), (k > 0), 
s*e°*9 = 9> s* e 1 = s s , s* £ k ( x ) = / ^ (fc_1) (x - u)s £ {u)du, (k>2). 
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2.2 Dynkin's formula for smooth subexponential functions 



Let us recall that taking expectations in the well-known Dynkin formula gives: 

Eg(X t )=g(0)+ [ E (Lg(X u )) du = g(0) + 1 ( E (Lg(X at )) da, (2.5) 



valid if g € . Iterating (|2.5I) . one obtains the following expansion for g € C?»+2 ( e .g^ see |24[ Proposition 9]): 

Eg(X t ) = 5 (0) + kl^ 9 ^ + J (1 " a) " E ( Ln+1 9( X -t)) da. (2.6) 



k=l 

Expansions of the form (|2.6p are then called iterated-type Dynkin expansions. For our purposes, it will be useful to 
extend (|2.6[) to subexponential functions. The proof of the following result can be found in Appendix El 

Lemma 2.1. Let v satisfy \2.I$ , and let g £ (j 2n + 2 fr e SU ch that 

limsupe-^ 1 !^?/)! < oo, (2.7) 

for any < i < 2n + 2. Then, \2. 6\) holds true. 
Below, let 

b := - / c e (u)v{du), bi :=b- / u(c e (u) - li^< 1 )iy(du), 

JR JR 

bi := cr 2 /2, 03 := — f u 2 c e {u)v{du), and 64 := / c e {u)v{du). 

2 Jr Jr 

Note that all these constants depend on e > 0, but this is not explicitly indicated for the ease of notation. In order to 
work with the iterated infinitesimal generator L k appearing in (|2.6|) . the forthcoming representation will turn out to 
be useful (see [13j Lemma 4.1] for its verification^) . 

Lemma 2.2. Let K, k = {k = (k , . . . , fc 4 ) g N 5 : k H h fc 4 = k} and for k e JC k , let £ k := k x + 2k 2 + 2k 3 . Then, 

for any k > 1 and e > 0, 

L k 9i*)= E IlftC£)B^g{x), (2.8) 

where 

!{ k 3 fe 4 \ 

/ .9 (M x + J2fo w 3 + J2 Ui dn K, > tih + ki> 0, 
5 ^(x), */fc 3 = fc4 = 0, 

and where the above integral is with respect to the probability measure 

rf7r k, e = H -£- c e( w 3) w ]v(dwj){l - /3j)df3j Y[ —c e (u i )i'(du l ), 

j=l 3 i=l 4 

on M. ks x [0, l] fc3 x R ki (under the standard conventions that 0/0 = 1 and Ili=i = 

Remark 2.3. The expansion (|2.8I) holds true for (possibly unbounded) functions g 6 £<2fc+2 sa tisfying (|2 . T[) for any 
< i < 2fc + 2. 



3 For convenience wc switched the role of 63 and 64. 
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3 Small-time expansions for the tail distributions and option prices 



In this section, we derive the small-time expansions for both the tail distribution V (Z t > z), z > 0, and for the call- 
option price function E (e z+Zt — l) , z ^ 0. With an approach similar to that in (T3J Theorem 3.2], the idea is to 
apply the following general moment expansion (easily obtained by conditioning on the number of jumps of the process 
XI introduced in Section |2~T|) : 



Ef(Z t ) = e-^Ef(U t + X[) + e- x ^ f] M^E/ ( C/ t + X t e + ^ ef ) (3.1) 

k=n+l ' \ i=l / 



e 



-x E t \ ~* 



^ k) 

fe=i 



-Ef[U t + XZ + J2Zi), (3-2) 



where {ff }i>i ar e the jumps of the process X s . We shall take f(u) = f z (u) := 1{ U > 2 } in order to obtain the expansion 
of the transition distribution and f{u) = f z (u) '■= (e z+u — 1)+ in order to obtain the expansion of the call-option price. 
To work out the terms in (|3.2[) . we use the iterated formula (|2.6p . while to estimate the terms in (|3.1[) , wc assume that 
the underlying stochastic volatility model U satisfies a short-time "large deviation principle" of the form: 

limtlnP([/ t > u) = -l-diu) 2 , (u>0), (3.3) 
t— >o 2 

where d(u) is a strictly positive rate function. In Section [3.41 we review conditions for ()3.3|) to hold. 

The expansions provided in the sequel will hold uniformly outside a neighborhood of the origin. Concretely, for 
fixed zo > and e > 0, the term O etZo (P) denotes a quantity, depending on z, e, and t, such that 

sup sup t^ 3 \O e . Z0 (t j )\ < oo, (3.4) 

0<t<t o |z|>z 

for some to > 0, small enough, depending itself on e and zq. 



3.1 Expansions for the tail distributions 



We first treat the case f z (u) := l{ u > z y The following expansion for the tail distributions of Z (whose proof can 
be found in Appendix [Bj holds true. 

Theorem 3.1. Let zq > 0, n > 1, and < £ < zo/(n + 1) A 1. Let the dynamics of Z be given by U.2\) and il.3]) . and 
the conditions \2.3\ )- §2J$ and \3. 3\) be satisfied. Then, there exists to > such that, for any z > zg and < t < to, 

V(Z t >z) = e~ x ' t J2^A^)- + O^ Z0 (t n+1 ), (3.5) 

3 = 1 J ' 

where 

k=l ^ ' 

with f k , z {y) ■= jTL y s* k (u)du. 

The expression p.5j) is not really satisfactory since the coefficients Aj t t are time-dependent and so the asymptotic 
behavior of the tail probability ¥(Z t > z) as t — > is unclear. In order to obtain an expansion of Aj it , we can further 
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obtain an iterated Dynkin expansion for Kg(Ut, Y t ) and a suitable moment function g. Indeed, assuming for simplicity 
that and arc independent, (U,Y) is a Markov process with infinitesimal generator 

„ / x , ,9o a 2 (y)d 2 q , »dq r ) 2 (y)d 2 q ,„ „, 

C 9 (u,y) = MJL + -f^ + «(y)^ + Xf^, (3.6) 

for g £ C 2 . Ito's formula and induction imply that 

n ifc +n+l rl 

E 5 (C/ t ,y t )= 5 ( Uo ,2/o) + X]T? £fe 5( u o,2/o) + — - / (l-a) n E{£ n+1 g(U aU Y at )}da, (3.7) 
fc=i ' n ' ■'° 

for any function g such that C k g(u, y) is well-defined and belongs to Cf, for any < k < 2n + 2. As in the case of 
the infinitesimal generator of X, one can view the operator (|3.6|) as the sum of four operators. However, given that 
in general those operators do not commute, it is not possible to write a simple closed-form expression for C k g(uo,yo) 
as is the case for X. Nevertheless, the following result gives a recursive method to get such an expression when 
fJ,(y) = —<J 2 {y)/2 and g(u,y) = h{u) as needed here . 

Proposition 3.2. Let the dynamics ofU andY be given by hi. 2)) and M.3\) with independent W 1 - 1 ^ andW^ 2 ' and with 
C°° deterministic functions a, a 2 , and 7 2 . For h,h £ C 2 (R), let 

C 1 h(u):=h"(u)-h'(u), C 2 h(y)~^-h"(y) + a(y)h'(y). (3.8) 
Then, for h £ (J 2n + 2 an d g(u,y) :— h{u), the infinitesimal generator h3. 6\) is such that 

k 

C k g(u,y) = Y J B k {y)C{h{u), for < k < n, 

3=0 

where B k (y) are defined iteratively as follows: 

B$(y) = l, B k (y) = 0, Vj g {1, . . . , k}, k > 1, 
B k (y)=£ 2 B k - 1 (y) + -^-B k ll(y), 1 < j < k, k>l. 

Proof. The proof is done by induction. □ 



Using the previous result, we can easily find conditions for the iterated formula (|3.7[) to hold. To this end, let us 
introduce the following class of functions: 

C[ l = { P eC n : \p {l) (x)\ < M„(l + \x\), for all < i < n, and 
for some Ai n <oo independent of a;}. 

tion [3 

Then, Ij3.7\ ) is satisfied for g(u, y) := h{u) whenever h £ C 2n ^ 



Corollary 3.3. In addition to the conditions of Proposition \ 3.2\ let 7 £ C ; and let a,a,j £ Cf, for any k > 0. 

" 1+2 

y 6 



Proof. Using Ito's formula and induction, (|3.7p will hold provided that 

'* dC n g(U s ,Y s ) 
du 

are true martingales. For this, it suffices that 



a(Y s )dWP, and f ^f'^ ^dW? 
Jo dy 



-q^-v] ds < 00, and E J (-q^-"/J ds < Vi > °- 
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Let us recall that since a and 7 belong to C ; °, 



supElY^™ < 00, 



(3.9) 



for any t > and m > 1 (this is similar to [261 Problem 5.3.15]). Hence, given the representation of Proposition 
it suffices to show that for some constants < 00 and non- negative integers r": 



(b?)«(iO <^(i + M) r? > 



(3.10) 



for any i, n > and < j < n. This claim can again be shown by induction since, given that it is satisfied for n — 1 
and using the iterative representation for in Proposition 13.21 



(5?) (<) (V)| <E 



«=o 



+ (a)W (B^f"^ + i (a 2 ) (£) (B^lf- l) 



which can be bounded by A4™(1 + |2/|) r '" since, by assumption, a 2 , a, and 7 2 belong to C, fe , for all k > 0. 



□ 



The previous result covers the Heston model, 



dU t 



1 



(2) 



(3.11) 



as well as the exponential Ornstein-Uhlenbeck (OU) model: 



dU t 



2Yt dt + e Yt dW^ 1] , dY t =x(0-Y t )dt + vdW t 



(2) 



(3.12) 



Both models are used in practice. 

Let us now use Corollary 13.31 to obtain a second order expansion for E/i([/ t ). Omitting, for the ease of notation, 
the evaluation of the functions Bj at yo, we can 



write 



where 



Eh(U t ) = h(0) + Bldh(0)t + (BfdhiO) + B 2 £jh(0)) t 2 + <3(/ 3 ), 



B\ = \a 2 , B 2 = ^ (a o% + K) 2 ) + a a a' Q , B\ = ±a 4 . 



(3.13) 



Above, we set ao = "'(2/0)1 °o = a ' (2/0)1 ancl = a " (2/0)1 with similar notation for the other functions. For the Heston 
model (|3.11|) . the coefficients in f|3.13[) are 



ol VO R 2 VO n 2 X{V-yo) 

b '-y> S2 "Y' Si ~ — 2 — 



Similarly, for the exponential OU model (|3.12l) . 



*i = Y' fi 2 2 = V' B? = e2vo ^ + x{e - yo) ) 



A general polynomial expansion of transition distributions is as follows. 
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Theorem 3.4. With the notations and the conditions of Theorem \3.1\ and Corollary \3.3\ 



P(Zt >z) = e-^ Y^Uz)- + Oe, zo (t n+1 ), (3.14) 



where 



E L J ar ) L9 Aj2 Br m(yo)^T (£,.) )(o), 



p+g+r— j x 7 \m— 

setting fo, z (y) = o^d where the summation is over all non-negative integers p,q,r. 

Proof. It is enough to plug the expansion (|3.7[) into the coefficients of the first summation in (|B-5|) (See the proof of 
Theorem l3.1l in Appendix B) and rearrange terms using Proposition l3.2l Note that the last integral in (|3.7[) is bounded 
for fk,z £ C^°(M) from the representation in Proposition 13.21 and the estimates (|3.9l) - (|3.10p . □ 

As a way to illustrate the expansions, note that in the case of constant volatility (a{y) = j(y) = 0), 



B k k (y) 

Hence, 



2 \ fe 

, Sj( s )sO, Vj^t, *>0. 



%<•>=- E_,C,l r )(|) ^(a(/,.))(o). 



Remark 3.5. (i) Clearly, the expansion Q3.14p leads to an expansion of the form: 



\Z t > z) =J2 & ^ z h + °^o(t n+1 )- (3-15) 



Indeed, by expanding e Ae * in (|3.14p . we have 

fc 

J 



a k {z)=Y J (%o^){-^) k - j - (3.16) 



(ii) The coefficients Sfc(z) in p,15[) are actually independent of e (for e small enough) since they can be defined 
iteratively as limits of P (X t > y) as follows: 



aUz) = lim -P (Z t > z) , = lim 

t->o t k\ • •• 



3.2 Expansions for the call option price 



1 f fe 1 #1 

|P(Z t >*)-$>(*)^| 



For z ^ and t > 0, let 

G t {z) :=E(e z+z < - l) (3.17) 
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where Z is the jump-diffusion process given by (|1.2I) and (|1.3[) . We proceed to derive the small-time expansion of Gt 
as t 4- 0. We first consider the out-of-the-money case z < from which one can easily derive the in-the-money case 
z > via put-call parity (see Corollary 13.91 below). Throughout this section, we set 

/(«) = /,(«) := (e 2+u -l)+, 

and we also assume the following uniform boundedness condition: there exists < M < oo, such that 

< o-(j/) < M. (3.18) 

Remark 3.6. Under the uniform boundedness condition (|3.18[) , it is easy to see that Ee cC/t < oo, for some c > 2. Then, 
a proof similar to that given in Corollary [231 using the representation of Proposition 13.21 shows that (|3.7[) is satisfied 
for g(u,y) := h(u), whenever h <G (j 2n+2 is a subexponential function satisfying (|2.7j) . 



The next theorem gives an expansion for the out-of-the-money call option prices in terms of the integro-differential 
operators L e and L\ defined in (|2.1|) and (|3.8|) (its proof is given in Appendix [Cjl . 

Theorem 3.7. Let zq < 0. n > 1. and < e < —zq/2(ti + 1) A 1. Le£ the dynamics of Z be given by U.2\) and U.3)) . 
and the conditions of both Theorem \3.1\ and Corollary \3.3\ as well as i3.18]) be satisfied. Then there exists to > such 
that, for any < t < to and z < zq, 

n p 

G t (z) = e ~^Y,b ] (z)- + O E , Z0 (t n+1 ), (3.19) 
i=i J ' 

where 

j 



K v,q,r j 

p+q+r=j v^'^' / \m=0 



b ^y-= E („„ r ) Lq e[z2 Br m (yo)£?(f P ,*))(o), (3.20) 

fo,z(y) = 0, and 



/*,,(!/) := / UV + «)^ fc ( W )du = / {e/+y +u - 1) + S^ fc («)d«. 



Remark 3.8. By expanding e Ae * in (|3.19p . one obtains (|1 .6[) with the following coefficients 

fc 



bk(z):=J2[fjb j ^)(-^) k - j - (3.21) 

j'=i 

As it was the case for the expansions of the tail probabilities (see Remark |3. 51 - (ii) above), the coefficients bk(z) can be 
defined iteratively as certain limits of Gt(z) = E (e z+Zt — 1J , as t — >• 0, and, therefore, they are independent of e (for 
e small enough). 

To deal with the in-the-money case z > 0, note that 

E (e z+Zt - 1) + = E{e z+Zt - 1) + (e z+Zt - 1)_ 
= e 2 - 1 + E(e 2+Zt -1)_. 

The expansion of E(e z+Z * — 1)_ with z > is similar to that of E (e z+Zt — l) + with z < 0. Therefore, we obtain the 
following result. 

Corollary 3.9. Let zq > 0, n > 1, and < e < zo/2(n + 1) A 1. Under conditions of Theorem \3.1\ there exists a 
to > swc/i that, for any < t < to, z > Zq, 

n y-m 

G t (z) = e z - 1 + e^ V 6 m (z)— + Oe, Zo (^ +1 ), (3-22) 
* — ' m! 

m— 1 
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where 
with 

9k Av) ■= [ {e y+z+u -l)_s?{u)du. 

JWl 

Remark 3.10. The results of this section provide expansions for the option price of out-of-the-money (OTM) and in-the- 
money (ITM) options in non-negative integer powers of the time to maturity. However, as stated in the introduction, 
at-the-money (ATM) option prices are known to be expandable as fractional powers, at least in the first leading term 
(see Proposition 5], [UJ Proposition 4.2], and [3T] ). A natural question is then to try to understand the reasons 
for such radically different asymptotic behaviors. The key assumption that allows us to obtain non- negative integer 
powers is the smoothness of the Levy density s outside the origin. Concretely, if s were discontinuous or were not 
differentiable at a fixed log-moneyness value z, we would expect that the expansions for the tail probability F(Z t > z) 
and for the OTM call option price E(e fe+Zf — 1) + would exhibit fractional powers. This fact was pointed out in [2"5] 
for the tail probabilities and a particular type of pure-Levy processes. The same is true for K(e k+Zt — 1) + , in view 
of its representation as the difference of two tail probabilities (see (14) in [14]). In view of this observation, it is now 
intuitively clear that ATM option prices typically exhibit fractional leading terms. Indeed, at-the-money option prices 
E(e Zt — 1)+ correspond to the log-moneyness z — and the Levy density s, is discontinuous at for an infinite-jump 
activity Levy process. 



3.3 Other payoff functions 



One of the advantages of our approach is that it can be applied to more general payoff functions. Concretely, 
consider a function of the form: 

fz(u) ■= tp(u)l {u > z} , 

where tp £ C£°. One can easily verify that, under the conditions of Theorem 13.41 and for z > 0, 

n p 

Ef z (Z t )=e- x ° t J2aj(z)^ + Oe, Z0 (t n+1 ), (3.23) 



where 



p+q+r= j V' y ' / \ m =0 / 

with fo,z(y) = 0, and fk,z(y) '■= J* R fz{y + u)sl k (u)du = Lp(y + u)s~l k (u)du. Indeed, from the proof of Theorem 13.1 
(which is the key to Theorem I3.4|) . the only step that requires some extra care is to justify that 



hAy) -=K / Ay + u)s* k (u)du 

is in C°° and that sup y \ Jk\{y)\ < oo. This is proved by verifying (via induction) that 

+ u)s* e k (u)du + \A EMr^- 1 -^)^- 1 ' * #(* - y). 

-v i=0 
Similarly, under the stronger conditions of Theorem 13.71 one can easily consider payoff functions of the form 

/*(«):= ¥>(«)!{«>-*}! (z<0), 
with <p £ C°° such that |(p"'(u)| < M^e u for some constant Mj < oo and all j > 0. 
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3.4 On the short-time large deviation principle for diffusions 



Large deviation results of the form (|3.3|) have recently been developed for different stochastic volatility (SV) models. 
For instance, for uncorrelatcd SV models, [17] shows (|3.3[) under the following conditions: 



• 


The function a is uniformly bounded and Lipschitz continuous in K. 


(3.24) 


• 


3 M 2 > Mi > 0, s.t. < Mi < a(y) A j(y) < a(y) V j(y) < M 2 < oo. 


(3.25) 


• 


cr, 7 € C°°, and o~(y) — > a±, 7(1/) —t 7±, as y — >• ±oo. 


(3.26) 


• 


cr and 7 are diffeomorphisms with cr > and 7' > 0. 


(3.27) 


• 


3 y c e R, such that a" > 0, 7" > for y < y c , a" < 0, 7" < 


(3.28) 




for y > y c and cr' V 7' < M < 00 for some M > 0. 




• 


7(er~ 1 (u)) 

The function u i— > is non-increasing. 


(3.29) 



We refer to [17] for an explicit expression for the rate function /, which is not relevant here. The Heston model p. lip 
(even with correlated Wiener processes and W^) was also considered in [TB] and j!8j . 

4 Particular examples and analysis of the parameter contributions 

In this section, we shed some light on the specific forms of the general expansions obtained in the previous section 
for some specific models. We also analyze the contribution, to the call option prices, of the various parameters of the 
model. 

4.1 Exponential Levy models 



Let us start by considering the exponential Levy model, which can be seen as a particular case of the jump-diffusion 
models (|1.2[) and <|1 .3|) by setting a(y) = 7(2/) = and a(y) = a. In this setting, the log-return process Z is a Levy 
process with generating triplet (b, cr 2 , v). Then, the following expansion for the out-of-the-money call option price holds 
true (see also [14] for an alternative derivation). Note that the process Z in this section is the same as the process X 
in Section 2 and therefore all the notations are transferred accordingly. 

Corollary 4.1. Let zq < 0, n > 1, and < e < — zo/2(n + 1) A 1. Let Z = (Zt)t>o be a Levy process with triplet 
(b,a 2 ,v) satisfying \2. S\) - [S^ . Then there exists to > such that, for any z < zq and < t < to, 

G t {z) = e-^j2c 3 {z)- ] +O e ^{t n+l ), (4.1) 

where 
with 

hkAv) ■= I {e z+u+u -l).sf{u)du. 

JR 

As in Corollary 13.91 we also have the following result for the in-the-money case. 
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Corollary 4.2. Let zq > 0, n > 1, and < e < Zo/(n + 1) A 1. Then, there exists to > such that, for any z > z$ 
and < t < to, 

" p 

G t (z) = e'-l + e-^ £ c,-(z) - + O £ , zo (t" +1 ), (4.2) 
j'=i J ' 

w/iere 

: = E Q^r^MCO), ftfc,.(l/) := jf (e W " - l)_^*(u)du. 
4.2 Second-order expansions 

Here are the first two coefficients of (|3.14[) for e > small enough: 

ai(z) = Bg(j/ )A,*(0) = / f z (u)s £ (u)du = / s(u)du; 

JR J z 

a 2 (z) = 2L e {f 1 , z ){0) + 2 Bl(y )C 1 (f 1 , z )(0) + / 2 ,*(0) 

p— l,q— 0,r— 1 j>— 2,q— 0,r— 

1 

2 , 




s'O - - /3)d/3u 2 s e (u)du 

<r 2 (yo)(s'(z) + s(z)) + / l{„ 1+U2 > 2 }S £ (iti)s E (it 2 )ditidit2- 
Jk 2 



The corresponding coefficients for (|3.19p are obtained as above with f z (x) = 1{ X > Z ] replaced by f z (y) '■= (e z+y — 1)+ 
with z < 0. Hence, for e > small enough, 



/ 

JR 2 



h(z)= / (e 2 +"-l) +S ( U )d M ; 

/oo 
e z+ "s(u)du 
-Z 

y+u 1+ u 2 _ y + s^ Ul )s e ( U2 )d Ul du 2 (4.3) 
(1 - p) ( [ e z - pu+w s(w)dw + s(~z - j3u)) d/3u 2 s £ (u)du. 

\J-z~l3u J 

In the previous expressions one can substitute c £ (y) and c £ {y) by lo< |y| <e an d l|y|>ei respectively 

Remark 4.3. Combining (|1.6[) . (JT77J) , and the expression for 6i(z) = b\(z) above gives the following expansion for the 
price function of the out-of-the-money call option near the expiration T: 

C(t,s) = Ke- r(T - t] G T -t{\ns - In if) 

= {T-t) [ (se u - K) + s{u)du + O eMs/k) ((T - t) 2 ) . (4.4) 



In particular, the leading term for out-of-the-money call option prices is "dominated" by the jump component of the 
model regardless of the underlying continuous stochastic volatility component U. 

Remark 4.4. Given that in (|3. 131) . Bq(yo) and B\(yo) depend only on o-(yo), it is interesting to note that the first two 
coefficients in our expansions (|3 . 14[) and (j3. 19[) coincide with the coefficients corresponding to an exponential Levy 
model with volatility a — cr(yo). In fact, the initial values of a and 7 begin to appear with the coefficients 03(2) and 
63(2) through the coefficients B\ in (|3 . 13|) . 
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4.3 Parameters contributions 



For simplicity, let Sq = 1 and r = 0. The first term on the right-hand side of (|4.3[) indicates the main contribution 
of the spot volatility cr(y ) in the call price. Indeed, recalling (|1.8p and the subsequent interpretation of — z as the 
log-moneyness k, one can say that the spot volatility has the effect of increasing the call price by 

a 2 (y )e K s(K)j (l + o(l)). (4.5) 

This observation was first pointed out in |14] for an exponential Levy model. 

As explained in Remark 14.41 the specific information present in the stochastic volatility model manifests itself 
up to an 0(t 3 ) term. Indeed, by analyzing each of the terms of the expansion for 63 given in (|3.20l) . the spot drift 
a(yo) and the diffusion component 7(2/0) of the underlying stochastic volatility factor Y are only felt through the term 
Bf(y )e K s(K)t 3 /2, which, in light of (I3T3)) . takes the form: 

H2/oK(yo) + W(yo)f) + a(j/oM2/oV (yo)) e K S (n) |. (4.6) 

For instance, in the Heston model (|3 . 1 1 1) . cr(y) — ^/y and the effect of the stochastic volatility starts to being felt 
through the third order term 

x{9 ~ ya) e* S (K)t\ (4.7) 

So, when the volatility is high (low), relative to its average value 6, a mean-reversion speed of x will decrease (increase) 
the call price by the amount (|4.7[) . This result is intuitive as the larger the speed is, the faster one would expect the 
volatility to come back to its average value when this is high, hence resulting in a discount. In a Heston model, the 
effect of the volatility of volatility v appears up to a fourth degree term of the form 

Aoe K (*"(«) -*'(«)) J V 
In contrast for the exponential OU model (|3.12|) . the quantity (|4.6[) simplifies to 

e^( X (0-y ) +V *) eKs{K)t3 

and both the speed x °f mean reversion and the volatility of volatility v already appear in the third order term. In 
spite of appearing in a third-order term, in some situations the contribution of the stochastic volatility model could be 
significant such as under a fast-mean reverting framework where the value of x is high. This phenomenon has already 
been suggested in the literature, partially motivated by empirical evidence (see, e.g., |19j). 

By analyzing the terms in (|3.20j) . one can see that in addition to the term (|4.6|) (which clearly involves the 
information of a at yo), the contribution of the volatility a(-) to the third-order term 63 also appears in the following 
term(s): 

(|^ 4 (yo) («"(«) + *'(«)) + - 2M'(«))) e"~. (4.8) 

Identifying the effect of each model component at each level of resolution t,t 2 , . . . , might be more enlightening 
than just explicitly writing down each of the terms. For instance, one can say that under a Heston model with 
independent Levy jumps, the first and second-order terms are the same as those of the underlying exponential Levy 
model with volatility o~(yo) = yo] while the third order term is the superposition of the corresponding third order term 
of the underlying exponential Levy model with volatility o~(yo) = t/o and the term (|4.7[) . From here, we can write the 
third order term of the underlying exponential Levy model using the expansion of Theorem 13.71 or using the following 
representation obtained in [TJ] (see the proof of Proposition 2.2 therein): 

e Xt - e K ) = ¥(X; >k)- e K P(X t > k). (4.9) 
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Here, X represents a Levy process with triplet (6, <j 2 {ya), v), while X* represents a Levy process with triplet (6*, a 2 (yo), v*), 
where 

1 



b=-^v 2 (yo)~ J (e x -l-xl^Mdx), v*{dx) = e x v{dx), (4.10) 
b* = \o 2 {w)~ I (e a -l-al w <i)i/(da;)+ f x (e x - 1) v{dx). (4.11) 

& J~oo J\x\<l 

Next, let dj(K;b,a,u) be the j -order term in the tail expansion of a Levy process X with triplet (6, a 2 ,v), i.e., the 
dj's are such that 

P(A t > k) =J2dj(n;b,a,iy)- +0(t n+1 ). 

3 = 1 J ' 



Then, the third-order term of the option price (|4.9[) takes the form 

/ \ t 3 

(d 3 (/e; b*,a(y ),v*) - e K d 3 (n; 6, a(y ), v)j ^ , 

which in turn allows us to write the third-order term of the model (|l.ip - (|1.3[) with the Heston specification (|3.1ip as 

d 3 (n; b*,a(y ), v*) - e K d 3 (n; b, a(y ), v) + M ° ~ Vo) e K S ( K f) |j. 

We can further disentangle the volatility effect from the second and third order terms by using ()4.5j) and (|4.8|) . 
Specifically, under the Heston model (where a(y) = y 7 ?/), the second-order term admits the representation 

(d 2 {n-b* ,0,0 - e K d 2 (K;b,0,is) + y e K s( K )) (4.12) 



while the third order term can be written as 

*! 

3! 



~( d 3 (K;b*,0,v*)-e K d 3 ( K ;b,0,v) 



+ ~y (sf(K) - 2b e s'(n) - 2 s (n)X e ) eA . 



(4.13) 



In these expressions, b* and b are computed as in (I4.10l) - (|4.11l) but setting a(yo) — 0. Note also that the term 
— s(n)e K \ e t 3 /2 in f|4. 1 3[) comes from expanding e~ Ae * in the second-order term e~ Xet b2(z)t 2 /2, and isolating the contri- 
bution of cr(j/o) in the latter term (see (|4.5I) V The expressions for d 3 (n; b* , 0, v*) and d 3 (n; b, 0, v) can easily be written 
from the expressions in |13| (see Theorem 4.3 therein) or from Theorem 13.41 above. 

Alternatively, one could also use (14.12[) and (|4.13l) to develop the second and third order SV corrections to call 
option prices of pure-jump Levy models. Concretely, the following expansions hold for the model (|l.l[) - (jl.3[) with the 
Heston SV specification (13. 11[) : 

E (St - e K ) + = E (e x * - e K ) + + a 2 (y )e K s( K ) t - + 0(t 3 ), 

= E(e x t~e K ) + +<j 2 (y )e« S ( K )^ 

+ ( 3x(6 y ^ eKs{K) + \ yl {s " {K) + S ' (K)) eK 

+ | Vo (sf («) - 2b £ s'(n) ~ 2 S ( K )X £ ) | + 0(t 4 ). 

Above, X denotes the pure-jump Levy model underlying S 1 , i.e., A is a Levy process with Levy triplet (6,0, v) 1 with b 
chosen as in \1A\ . The call option prices E (e Xt — eK ) + can be computed using Fourier inversion formulas or, in some 
particular cases, via analytic closed form formulas. 
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5 Asymptotics of the implied volatility 



Using the leading term of the time-i price for the out-of-the-money call option as computed in the previous section, we 
now obtain the asymptotic behavior of the implied volatility a{t\ s) near T. This is defined implicitly by the equation 

C{t,s) = C BS (t,s;a(t;s),r), (5.1) 

where OBs(t, s; a, r) is the classical timc-i Black-Scholes call-option price corresponding to an interest rate r, a volatility 
<7, and time-i spot price s. First, we need the following well-known result (see, e.g., Lemma 2.5 in [2"U]). 

Lemma 5.1. Let Cssit, s; a, r) be the classical Black-Scholes call price function. Then, as t j' T, 

1 Ka\T-t)^ ( {\nK-\nsf \ 
° BS ^ 55 T) = 7^ (In A -Ins) 2 ° XP ar»( T -t) J (5 " 2) 

/ InA-lns r(lnA-lns)\ 

x exp v 2 — + V 2 j +i H*' s ; CT ' r )- 

The remainder satisfies 

\R(t, s; a, r)| < M(T - tf' 2 exp (- ^f" jffi ) , (5.3) 

where M = M(s,o~, r, if) is uniformly bounded if all the indicated parameters vary in a bounded region. 

The next result gives the asymptotic behavior of a(t,s). This has already been obtained for a pure-Levy processes 
(see, e.g., and [H]) and is presented here for the sake of completeness. 



Proposition 5.2. Let the dynamics of Z be given by (QUI) and S1.3\) . and the conditions of both Theorem \3.1\ and 
Corollary \3.3\ as well as 13.18]) be satisfied. Let &(t;s) be the implied volatility when the time-t stock price St is s. 
Then, as t j" T , 

, N (InA-lns) 2 , N 

^ S ^ -2(T-t)HT-ty ^ 

Proof. Let k := In A/s be the log-moneyness and let r := T — t be the time-to-maturity, where, for simplicity, we write 
a{t) instead of a(t;s). Using the leading terms in (|4.4I) and (|5.2[) , we obtain that as t j T: 

tu{s,K)~ v{s,K)a 3 (t)T 3 / 2 exp ( -— + , (5.5) 



where 



2a 2 {t)r a 2 {t) 



f 1 K 

u(s,K)= / (se u — K ) + s(u)du, v(s, A) = —= — e =, 

JR V27T K 

.1/2 



Assume that limsup t ^ T <j(t)r ' = c G (0, +oo), then limsup t -|. T a(t) = +oo and, thus, 

hmsup (a(t)r^f exp f-^- + ^) = c 3 exp (-^ ) 



So the right hand side of (|5.5p does not converge to while the left hand side does, which is clearly a contradiction. 
Now if limsup t ^ T a{t)T X / 2 = +oo, then limsup t ^ T a(t) = +oo and, thus, 



limsup ((j(t)r ' ) exp 



2a 2 {t)r & 2 {t) 



Again we obtain the same contradiction. Therefore, we have limsup t ^ T (r(t)r 1 / 2 = 0. Then, (|5.5[) can now be equiva- 
lently written as 



u(s,K) 

2a 2 (t)T ' & 2 (t) ' "'"V" v "'' J ~' y 'w(s,A)' 



exp(- " + _g_ + 31n (a(t)T 1 / 2 ) -lnr 
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Hence, as 1 1 T, 



lim (- tt4u^ + + 3 ln f °{ty /2 ) - In t - In 4^ ) = 

tfT V 2er 2 (t)r cr 2 (i) V w / v(s,K)/ 



Finally, since we obtain (|5.4[) . □ 

Remark 5.3. 

1. As seen in Proposition 15.21 the leading order term of the implied volatility near expiration is "model free", i.e., 
it does not depend on any of the model parameters. 

2. As discussed in Remark POI the second order expansion of OTM call option prices coincides with that of a purely 
exponential Levy model with volatility parameter a equal to the spot volatility <j(Y t ). Thus, the second-order 
expansion for the implied volatility given in |14) applies: 

a 2 (t;s) = v (t-k) (l + v^t; k) + o (j^xj) - ( r -> 0). 

where re and r denote respectively the spot log-moneyness re := log K / s and time-to-maturity r :— T — t, while 
vq and Vi are given by 

ire 2 

v (t; re) — 



-rlogr 

vi{t;k) = bgW iog ^V' 2 1 (e " - eK )+ s H^ io g 3/2 (^)) • 

3. In the very recent manuscript |22j . the authors give a blueprint to generate expansions of arbitrary order for the 
implied volatility i7 2 (t; s). Interestingly, such expansions are determined exclusively by the leading order term of 
the option price expansion, meaning that the stochastic volatility correction is not relevant for implied volatility. 

6 Small-time expansions for the Levy transition densities 

In this part, we revisit the important problem of finding small-time expansions for the transition densities of Levy 
processes. This problem has recently been considered in [37] and also in [TS]. As in Section |2~T1 we consider a general 
Levy process X with Levy triplet (6, tr 2 ,^). It is well-known that under general conditions (see, e.g, [IS] and [3~4"]): 

hmlf t (x) = s(x), (i,t0), (6.1) 

where ft denotes the probability density of X t and s is the Levy density of v (both densities are assumed to exist). In 
many applications, the following uniform convergence result is more desirable 



lim sup 

\x\>T) 



\fti X ) - S ( X ) 



0, (6.2) 



for a fixed 77 > 0. The limit (|6.2[) is related to the general expansions of the transition densities: 

Mx) = f2^P-t k + t^O v (l), (6.3) 
fe=i 



4 We thank an anonymous referee for bringing our attention to this manuscript. 
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which is valid for any \x\ > 77 and < t < to, with to possibly depending on the given 77 > and n > 0. Above, and as 
in (|3.4[) . Ot;(1) denotes a function of x and t such that 

sup sup I < 00. 

0<t<* |x|>!7 

Note that (|6.2[) follows from (|6.3p when n = 1 and afc(x) = s(x). 

The expansion ([6.20 was first proposed in [37] building on results of [28] , where the pointwise convergence in (|6.1|) 



was obtained. In both papers, the standing assumptions on the Levy density s of the Levy process X are: 

z 2 s(z)dz > 0; (6.4) 



-v 



(ii) S eC°°(R\{0}); (6.5) 

r is'Mi 2 

(iii) / 1 ) < 00, for all 77 > 0; (6.6) 

J\z\> v S K Z ) 

(iv) There exists h e C°° such that fc(z) = 0(z 2 ) (z -> 0), (6.7) 



h(z) > if s(z) > 0, and 



^-h(z)s(z) 
dz 



2 1 



-dz < 00. 

s(z) 



Condition (j6.4p is used to conclude the existence of a C°° transition density f t (see [35] Chapter 5]), while (|6.5[) - (|6.7[) 
are needed to establish an estimate for the transition density using Malliavin Calculus. However, the method of proof 
of [37) is not convincing and can only yield the first order expansion in (|6.3p (see the introduction of |13j for more 
details). Recently, [T3] obtained (|6.3[) under the following two hypotheses: 

j v , k := sup \s (k) (x)\ < 00, (6.8) 

\x\>Tj 

limsup sup \ft (%)\ < o°, for all k > and for all 77 > 0. (6-9) 

t\0 \x\>rj 

Condition (|6.8p is quite mild but condition (|6.9[) could be hard to prove in general since closed-form expressions of the 
densities ft are lacking. Nevertheless [13] shows that condition (|6.9[) is satisfied by, e.g., the CGMY model of [9] or 
Koponen [27] and by other types of tempered stable Levy processes (as defined in [35]). 

In this section, we show that (|6.9p is not needed to obtain (|6.3p . See Appendix [C] for the proof of the following 
result. 



Theorem 6.1. Let rj > and n > 1, and let the conditions |fi.^| /'-( (7Pj) be satisfied. Then, i6.3\) holds true for all 
< t < 1 and \x\ > 77. Moreover, there exists £0(77,72) > smc/i t/iat /or aZ/ < e < £o ; i/ie coefficients au admit the 
following representation (which is moreover independent of e for any < e < £0): 



ak( 



j=l ~^ J/ i=i 



(*) : =E \ X]l ! X ■')'• ( 6 - 10 ) 



where Si. x (u) := s* l (x — it). 

Remark 6.2. Combining the proofs of Theorem 13.11 and Theorem 16. 1[ it is possible to obtain a small-time expansion 



for the jump-diffusion model (|1.2[) - (|1.3| l assuming, for instance, that the stochastic volatility model admits a density 
function dt satisfying the small-time estimate: 

sup dt{x) < M PiV t p , 

\x\~>T) 

for any p > 1 and < t < to(p, rf) and some constant M p ^ < 00. 

Acknowledgments: The authors are grateful to two anonymous referees for their constructive and insightful com- 
ments that greatly helped to improve the paper. 
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A Proof of Lemma 12.11 

Let us show (|2.6[) for n = 1 (the other cases are easily obtained by induction). First, applying Ito's lemma ([ 
Theorem 1.4.56]), 

g{X t )=g{0)+ [ Lg(X u )du + a [ g'(X u )dW u 



(g(X u - + z)- g{X u ~))^j,(du,dz), 



where L is given by (| 1 . 10|) . One can easily check that Lg{x) is well-defined in light of the continuity of g^ and (|2.7I 
Indeed, there exist constants Mj,i = 0,1,2, such that |</ J )(:r)| < MjesM, for all x, and thus, 



g(x + z)v{dz) 



\z\>\ 



< M I e^v{dz)e^ x \ 

\z\>l 



(g(x + z) - g{x) - g'(x)z)v{dz) 



z\<l 



< M 2 e* 



kl<i 



l v(dz)e^ x K 



Next, we show that the last two terms of the expansion of g(X t ) above are true martingales. Indeed, it suffices that 



E / \g'(X u )\ z du<oo, 
Jo 

E / \g{X u + z) - g{X u )\v{dz)du < oo, 

JO J\z\>l 

e[ [ \g(X u + z)-g{X u )\ 2 v(dz)du<oc. 

Jo J\z\<l 



(A-l) 
(A-2) 

(A-3) 



Using (|2.7|l and the continuity of g', there exists a constant M > such that 

E f \g'(X u )f du < M [ Ee clXul du<M [ Ee cX "du+ [ Ee- cX "du < oo, 
Jo Jo Jo Jo 

for any t > 0. Similarly, setting B = {z : \z\ > 1}, (|A-2I) is satisfied since 



g'(X u + w)dw 



>o 



v(dz)du 



Ef f \g{X u + z)-g(X u )\v{dz)du< 
Jo Jb Jo Jb 

<M ( Ee clXul du I I e cw dwv(dz) < oo. 
Jo Jb Jo 



k-l 



Also, setting B = {z : \z\ < 1}, 



f f \g{X u + z)-g{X u )\ 2 V {dz)du<E f f f \g' \X U + zp)\ 2 dpz 2 v(dz)du 
lo Jb Jo JbJo 

< f Ee c ^du [ f e cW d(5 z 2 iy(dz) < oo. 
Jo Jb Jo 



We then have that 



Eg(X t ) = 5 (0) + E / Lg(X u )du, 



which leads to (|2.5[) . provided f*E \Lg(X u )\ du < oo. The latter is proved using (|2.7[) and arguments as above. 
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In order to obtain (12.61) for n = 1 by iterating (|2.5|) . we need to show that for any C 4 function g satisfying (|2.7I 



limsup e-^\{Lg) (i Hy)\ < oo, (A-4) 
for i = 0, 1, 2. To this end, we first note that 

(Lff)W(y) = bg^\y) + ^ +2 \y) + f (g^(y + z) - g^(y) - zg^^l^Mdz) 
for i — 0,1,2. Hence, it is sufficient to show (|A-4[) when i = 0. But, 

e-« |y| |L 3 (y)| < be-^\g'(y)\ + ^ e -^\g"(y)\ (A-5) 

+ e-fl»l / |«7(y + 2)-«7(w)k(<te) (A-6) 
J| z |>i 



e -Sly| 



+ «) - .9(2/) - (A-7) 



The limits of each term of the right-hand terms in (|A-5[) are trivially finite as \y\ — > 00 by the assumption (|2.7p . For 
the term in (|A-6[) . again by the assumption (|2.7|) and the continuity of gW, there exists M > such that, 

| ff W(y)|<MeSW i = 0,1,2. 

It then follows that 

e S I ^ I / \g(y + z) — g{y)\u{dz) = e — § Ivl / / + w)|eky i/(dz) 
J\z\>l J\z\>i Jo 



< M f (I e% w dw)v{dz) 



|z|>l V0 

= Ml e^v(dz) < 00, 
./|z|>i 

which immediately implies that 

limsup e"^ |y| / \g(y + z) - g(y)\v(dz) < 00. (A-8) 
lv|-n» J\z\>l 

Similarly, we can show that the limit as \y\ — > 00 of (|A-7j) is finite. Therefore, we can iterate (|2.5j) to obtain (|2.6j) for 
n = 1. □ 

B Proof of Theorem [3J] 

We analyze each term on the right-hand side of the expansion of E/(Z t ) given in (|3.1[l - (|3.2[) : 
(1) For any z > zq, we have 

Ef z (U t + XI) = F(U t + XI >z)< P(U t > z/2) + P(A t £ > z/2). (B-l) 
By our assumption (|3 . 3[) . there exists to(zo) > such that for any < t < to. z > zq > 0, 

^{Ut > z/2) < ¥(U t > z Q /2) < exp f- ^fH , (B-2) 
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which can be seen to be O zo (t n+1 ). Also, the second term on the right-hand-side of (|B-1[) is £yZ „(t n+1 ) in light of 
(|2.2|) by taking a := (n + l)/z and since < e < Zo/(n + 1) A 1. 



(2) The second term in (|3 . 1 1) is also O E:Z0 (t n+1 ) since / < 1 and clearly 

oo 

e" Ae * E ( X et) k /k\ ^ ( X ^) n+1 = 0(t n+1 ). 

fc=n+l 

(3) We proceed to deal with the terms in (I3.2|) . By the independence of U and X, 

E/ z lUt + Xt + J = EA,, (C/ t + A7) = E/ M , t (X t e 



; = 1 



where 



fk,M~(*e) «T(u)d« and / M , t (y):=E/ M (tf, + y) 



In particular, by the assumption (|2 . 3[) . 



sup 



< A^H^IU < A- 1 max %e/2 := T e < oo. 



It follows that /jb } a,t G C£°(R) and moreover, 



= E /S + tf) > and sup /g )t (tf) 



< r e , for any j > 0. 



We can thus apply the iterated formula (|2.6p to get 



n—k 



^kzAxi) = E 7f i*A I ,,t(o) 



i=0 



(n - fc)! y 



(1 - a) n - k E{Ur k+1 fk, z ,t( X£ a t)}da. 



It follows from the representation in Lemma 12.21 and (|B-4[) that 

sup f (l-a)"- i E(ir t+1 / iAt K t ))dQ<(X), 

z JO 

and thus the second term on the right hand side of (|B-5|) is O e . Zfl (t n ~ k+1 ). 



(4) Combining ([33]) . (pHj) and (|B^5]) . we obtain 



fc=i 



fc=l ' i=0 



k n Z^ J4, 



n+l\ 



i=i fc=i 

Using again the representation in Lemma 12.21 and (|B-4[) , it follows that 

Li- k f k ,zM = L e \®fk,z(U t + ■)] (x) = \- k L k £ 1 

and (13.51) is obtained. 



(B-3) 



(B-4) 



(B-5) 



□ 
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C Proof of Theorem 13.7 



We analyze each term in (|3.1j) and (|3.2j) through the following steps: 
(1) For z < 2 < 0, 



lf z {U t + XI) = E [e' +u * +x ' - lj + < E (e^+ x * l {c/t+Xf >_, } ) (C-l) 

1/2 



< (Ee 2l/t+2X ' P(t/ t + X t £ > 

< (Ee 2Ut Ee 2X ^ 1/2 (p(U t > -z/2)+P(X' > -z/2)^ P ', 
= e 1 ^' 2 (Ee 2Ut ) 1/2 (p(U t > -z/2) + nX! > -2/2)) ^ 

where "0 is the characteristic exponent of X e . Since M t := e 17 ' satisfies the SDE dM t = M t a(Y t )dW t {1) , and from the 
Burkholder-Davis-Gundy inequality, 

Ee 2Ut = E ^1 + J M s a{Y s )dW^A 

< 2 + 2E e u 'a(Y s )dW^ < 2 + 27\/ 2 E J e 2U °ds. 

By Gronwall's inequality, 

Ee 2U < <2e 2MH = O e , Zo (l). 
Therefore, the right-hand-side of (jCTT]) is of order o £ . Zo (t n+1 ) by (EH) and (pH]) . 

(2) The second summation in (|3.ip is also O £i 2 (i™ +1 ) since for any k > n+ 1, 

E/ z (C/t + X t e + V&) < e z Ee c/t Ee^(Ee« 1 ) fe < A7 fc e** (1) ( / e x s e {x)dx) . 

(3) To study the summation in (|3.2j) . recall that by the independence of [/ and X, for any 1 < k < n, 

where 

A,,,t(x) = E/ M (C/ t + ar) and 7 M (ar) = E/„ fx + ^ &j . 
Let us show that fk.z is C°° . Indeed, since 

f />OC 

/ M (x) - A £ ~ fc / / u t _ A ^ Ul ) dui TT s £ ( Ui )du e , 

JS.h-1 J- V* „u t -z-x V 7 7 o 
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and s £ G C£°, we have that 

/ e» +a+ ^«'a e («i) t i«iTT* e («/)d«/ 1 

k 



. / k \ k 

+ \J k / s £ I — }ut - z - x TT s e (ui)du£. 



Using induction, we see that 

/■ poo k 

fi%) = K k e z+x+ ^^s s (ui)d Ul l[s £ (ue)du e (C-2) 

JR J- JZe = 2 ui-z-x g_ 2 



i-2 / k \ k 



3=0 " m \ £=2 / fc2 

In view of (|2.3I) . there exists a constant M i£ < oo such that, for any i > 1, 

A' 



i-2 - k 

+ M h£ Xj k ^ / ^ s £ (u e )du e max 7^ 



(C-3) 



e/2- 

The right-hand side of (|C-3|) is integrable since Ee* 7 * = 1. By dominated convergence, we conclude that fk,z,t & C°°(l 
and also, 



fkl,ti x ) = E /£i(&t + x) , for all i > 0, and limsupe - *!*! f^ t (x) 



< 00, 



since c > 2. Thus, applying (|2.6j) gives 



4-1 j.n—K+1 rl 
m,z,tm = E - L ekzM + 7 TTj / (1 ~ «)"- t E{L»- t+1 A A ^ot)}*- (C-4) 

l=0 z - l n — K )- Jo 



To show that the last integral in (|C-4j) is bounded, we apply Lemma \2 .21 to get that 



'n — k + 1 



E((LrW/w)K t ))= E II 6 H t )E(B k ../M,*TO) 

ke/C„_ fc+ i i=0 



Thus, it is sufficient to show the boundedness of E£? k e fk,z,t(X^ t )-, for any 1 < < n and k = (kg, ■ ■ ■ , /C4) € /C„_fc+i. 
Indeed, noting that (|2.4p implies that 
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we have, for any x £ K and some constants K\, K2 < 00 



[0,l] fc 3 xI*3+*4 



J k.z.t 



< 



3 



[0,l] fc 3 Xl*3+ t 4 



Jk,z 



+ E fa™* + E u * 
j=i i=i y 

feg ^4 

U t + X + ^2 P] W 3 + E M » d7T k 
3=1 i=l 
it 



• .\/A ' Cr' / f e z+x+ ^*=i ue s e (ui)duif[s s {ue)dut 
Jm. k - 1 Jr e=2 

£ k -2 . k 

+ M ijS \~ k / ~\s E (u e )due max j js/2 
~k ii'-i ft 0<j<i 



j=0 " * 1=2 

= M x e x + M 2 < 00, 

where the third inequality follows from (fC^3]> . It follows that EB K J kjZit (Xg t ) is O 6 , Z0 (l), and so is EL^- fe+1 / M , t (X= t ). 
Therefore, the last integral in (|C-4|) is indeed O £j2o (i n_fe+1 ). 



(4) Plugging ()C-4j) into (|3.ip and (|3.2j) . and rearranging terms lead to 

Ef z (Z t ) = e -^Y,-J^k*A X t) + Os, Z0 (t n+1 ) 



fe=i 



e " Aet E S E (l ) ^r7M,t(0) + Oe^„(t" +1 ). 



7! 

3 = 1 J fe=l 

It remains to expand the coefficients 

£r fc /M,t(0) = L{- k [Ef k , z (U t + •)] (0) = \- h U- h [EfkAUt + •)] (0). 
Using the expansion (|3.7[) and Remark 13.61 we have 



(C-5) 



(G-6) 



fn-i+l rl 



(n-i+l)!J 

E|E^'(yo)4/MW 

i=0 ' Z=0 

tfi— 3+1 /■! 



da 



i=0 
n-3 i 



(C-7) 



(n-i + l)!j 

Finally, by combining (|C-5|) . (|C-6|) and (|C-7[) . it follows that 



(1 - a) n "% £ n - J ' +1 /fc,*(^t + x) da 



3 = 1 J fc=l 
j-n-j + 1 r 1 



n-3 



E^r" E b 'W £ '/m (o) 



i=0 



J=0 



+ 7 -rry / {l-a) n -m{Ll- k p^f k , z {U at + 

(n— 7 + 1 ! L L 



r/n 



+ O e ,z (^ +1 )- 



(C-8) 



(n-j + l)iy 

Finally, since the integral in (|C-8[) is O e . Zo (l), as seen from the uniform boundedness condition (13.18[) and the estimate 
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(|C-3[) . we obtain that 



EMZ t ) = e-" E^E (l) E I E 3(«>)4&.. ) (0) + O e> * (t" +1 ) 

3=1 -/ ' fc=l V / i=0 ' \l=0 / 



eS e Li j^(E s -^) £ ^](°)+ o ^o(* n+i )- 



7! * — ' \r>.q,r, 



□ 



D Proof of Theorem 16.11 

We only consider the case x > (the case x < can be done similarly by considering P(Xt < x)). Again, we start 
with the expression 



P(X t > x)= f Kt f{X e t > aQ +e-*'* £ ^f~ P ( X t + E & ^ x ) P" 1 ) 

fc=n+l ' V i=l / 



A- 



Bt(x) 



-A.t ' 



C t (x) 



fc! 

fe=l \ 1=1 



D t (x) 

Let / t e denote the density of Xf, whose existence follows from (|6.4[) . Given that 

k 



inxi + ±t i >x)='r t *,? {x i 



and that 

I £ e . —*k i „_\ i ^- i /■ i ^ _ 

e/2,0 A e 



sup |/ t £ * (x)\ < sup < ls/2 ,X k -\ 



one can interchange derivative and summation in (|D-1[) to show that for each t > 0, Ct admits a density ct, with 
moreover, 

00 j-k 00 / \ j-\k 

S np\c t (x)\=snpe-^ £ ^j? * (a;) < e "^^£ £ < A? 7e /2,o* n+1 - (D-3) 

fc— n+1 fc— n+1 

Also, in view of Proposition III. 2 in [28], there exists a real £o(?7, n) > such that for all < e < 6q and t < 1, 

sup < M(r7, £ )t"+\ (D-4) 

| a: | > rj 

where M (77, e) is some constant depending only on 77 and e. The last step of the proof is to deal with the terms in D t . 
Recall that 



i=l 



and 



zl fkAv) = x-^-iy^K^ * st'Hx y), 
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with 



( fc \ roo 

y + J2^^ x = A « fc / ■ 



Then, applying the iterated formula (|2.6[) . we get 

Ef k>x (X?) = ^ ^/m(0) + (^yy jf 



i=0 



s* E k (u)du. 



(1 - a)™-*E L™ +1 ~ k f k>x (X e at ) da 



Using the representation of L e in Lemma 12 .21 one can easily verify that 

d 



d.r 



sup 



< M n , fc , eo max n { 7s/2lfe }, 



(D-5) 

(D-6) 
(D-7) 



for some constants M n ^ k _ e < oo. Hence, one can pass d/dx through the integral and the expectation in the last term 
of (|T>5|) to get 

^-Ef k . x (X!) = -(\e)- k Y< ^ L %,M +t n+1 - k O E , k , n (l) 1 (D-8) 
x i=0 l ' 

where O e , k , n (l) indicates that sup,,, |O e ,fc, ra (l)| is bounded by a constant depending only on e, k, and n. Differentiating 
F(X t > x) in (|D-ip and plugging in (|D 3|) . (|D-4|) . (|D-8|I gives for any < £ < e and i < 1, 

n k n — k j 

f t (x) = e'^ E fc? E Jj L %A0) + * n+1 Oe,,(l), 

fe=l ' i=0 

where O e!?? (l) is such that sup t<1 sup^i^ |0 E)I7 (1)| < oo. Rearranging the terms above leads to 

f t {x) = e- x * 1 E c p {x)- + t n+1 O e , v (l), 



with 



The expression in (|6.10[) follows from the Taylor expansion of e~ Xct , using also that sup x |c p (a;)| < oo (a fact which 
follows from (|D-6|0 . Finally, the "constant property" of (|6.10[) . for any < e < Eq, follows from inversion. Indeed, 
given that a posterior 



k=l 



(D-9) 



holds true for any t < 1 and < e < eq, a k {x) can be recovered from ft(x) (independently of e) by the recursive 
formulas: 

ai (x) = \imjMx), 0fc (z) = lim^/ t (s)-£^t^, 2 < k < n. 



□ 
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